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We report the first observation of e + e~ annihilations into states of positive C-parity, namely p°p° 
and 0p°. The two states are observed in the tt + tt~-k + -k~ and K + K~ir + ir~ final states, respectively, 
in a data sample of 225 fb -1 collected by the BABAR experiment at the PEP-1I e + e~ storage 
rings at energies near \fs = f0.58 GeV. The distributions of cos^*, where 6* is the center-of- 
mass polar angle of the <j> meson or the forward p° meson, suggest production by two-virtual- 
photon annihilation. We measure cross sections within the range |cos(9*| < 0.8 of a(e + e~ — > 
p°p°) = 20.7±0.7(stat)±2.7(syst) fb and a(e + e"^ (jip ) = 5.7 ± 0.5(stat) ± 0.8(syst) fb. 

PACS numbers: 12.20.-m, 12.38.Aw, 12.40.Vv 



The process e + e~ — ► hadrons at center-of-mass (CM) 
energy y/s far below the Z° mass is dominated by annihi- 
lation via a single virtual photon with charge-conjugation 
parity C = — 1. The high luminosity of the B factories 
provides an opportunity to explore rare, low multiplic- 
ity final states with C = +1 such as those produced 
in the two-virtual-photon annihilation (TVPA) process 
depicted in Fig. 2] The TVPA process has been ig- 
nored in the interpretation of the total hadronic cross 
section in e + e~ annihilations as input to calculations Q 
of the muon g—2 and the running QED coupling a. We 
report the first observation of the exclusive reactions 
e + e~^ p°p° and e + e~^ <j)p°, in which the final states 
are even under charge conjugation, and therefore cannot 
be produced via single-photon annihilation. 

This analysis uses a 205 fb" 1 data sample of e + e~ 
collisions collected on the T(4S) resonance and 20 fb -1 
collected 40MeV below with the BABAR detector at 
the SLAC PEP-II asymmetric-energy B factory. The 
BABAR detector is described in detail elsewhere 0. 
Charged-particle momenta and energy loss are measured 
in the tracking system which consists of a silicon vertex 
tracker (SVT) and a drift chamber (DCH). Electrons and 
photons are detected in a Csf(Tl) calorimeter (EMC). 
An internally reflecting ring-imaging Cherenkov detector 
(DIRC) provides charged particle identification (PfD). 
An instrumented magnetic flux return (IFR) provides 
identification of muons. Kaon and pion candidates are 
identified using likelihoods of particle hypotheses calcu- 
lated from the specific ionization in the DCH and SVT, 
and the Cherenkov angle measured in the DfRC. Elec- 
trons are identified by the ratio of the energy deposited 
in the EMC to the momentum and by the shower shape; 
muons arc identified by the depth of penetration into the 
IFR. 

Events with four well-reconstructed charged tracks and 
a total charge of zero are selected. Charged tracks are 
required to have at least 12 DCH hits and a polar angle 
in the range 0.41 < < 2.54 radians. The momenta 
of kaon and pion candidates are required to be greater 
than 800 and 600 MeV/c, respectively. Among the four 
selected tracks, two oppositely charged tracks must be 
identified as pions, and the other pair must be identified 




FIG. 1: Two-virtual-photon annihilation diagram. 



as two pions or two kaons. Events in which one or more 
pion candidates are identified as an electron or muon are 
rejected (lepton veto). We fit the four tracks to a common 
vertex, and require the \ 2 probability to exceed 0.1%. We 
accept events with reconstructed invariant mass within 
170 MeV/c 2 of the nominal CM energy (Fig. 01. 

In the process e + e~ — > tt + tt~tt + tt~ , there are two pos- 
sible pairings of 7r + mesons with 7r~ mesons. However, 
only one combination appears in the kinematic region of 
interest (m w + w - < 2 GeV/c 2 ) for both pairs. We label 
the pion pair with CM momentum vector pointing into 
the hemisphere defined by the e~ beam direction tt + ttJ 
and the other as tt+tt^. Figure E^a) shows the scatter 
plot of the invariant masses of ft ttJ and tt+tt^ from 
e + e~ — * 7r + 7T~7r + 7r~ events, and Fig.^b) the plot of in- 
variant masses of K + K~ and 7r + 7r~ pairs from e + e~ — > 
K + K~7r + ir~ events. We observe correlations of masses 
in Fig.OJa) indicating the production of p°p° final states, 
and in Fig. I^b) indicating the production of 4>p° final 
states. 

To extract the number of e + e~ — > p° p° and <j>p signal 
events, we perform a binned maximum-likelihood fit for 
nine rectangular regions (tiles) in the two-dimensional 
mass distributions, as shown in Fig. [21 The signal box 
is the central tile (tile 5), defined by the mass ranges 
0.5 < m„+„- < 1.1 GeV/c 2 and 1.008 < m K+K - < 
1.035 GeV/c 2 . For e+e" -> K+K-jt+tt-, the expected 
number of events, n*, for each tile i can be expressed as: 
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FIG. 2: Distributions of the invariant mass (T(4S) data) for 
the a) 7T + 7r _ 7r + 7r _ and b) K + K~ ir + iv~ final states. The ac- 
cepted signal regions are indicated by the dashed lines. 
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FIG. 3: Scatter plots of the invariant masses of the two oppo- 
sitely charged pairs in the a) -k + it~-k + -k~ and b) K + K~-k + tt~ 
final states. The dashed lines indicate K + K~ /it + -k~ thresh- 
olds. The solid lines show the nine tiles used in the fit. 



m = ffS + ffN* + ff N p0 +f*B, (1) 

where S is the number of c/)p a signal events, is the 
number of 4>X background events, N p o is the number of 
p°X background events, and B is the number of resid- 
ual background events, in all nine tiles. The parameter 
ff is the fraction of events of type T that contributes 
to tile i. The signal fractions ff are modeled by Monte 

Carlo (MC) simulation and ff and /f are obtained 
from the cj>X and p°X background shapes, which are es- 
timated by fitting the projections of m K + K - and m^,- 
as described later. The residual background fractions ff 
are modeled by a linear function that can be expressed 
as 

fB _ AziAj/ifl + s p o{x l - x 5 ) + s^{yj - y 5 )} 
Ej=i AaSjAj/j 

where Ax* and Ay^ are the kinematically accessible di- 
mensions of tile i, Xi and yi are at the center of tile i, and 
SpO and s,/, are slopes obtained from the fits. A similar 
expression is used for the tt + tt^tt + tt~ case, where </> and 
p° are replaced with p° and pg. 

The background fractions are obtained by mass pro- 
jection fits which are confined to the central horizon- 
tal or central vertical <j> or p° resonance band. The ef- 
fect of neglecting the resonance width outside the central 
band, checked by smearing the background fractions in 
the central band into the adjacent tiles using the reso- 
nance widths obtained from MC, is found to be negligi- 
ble. The mass projections in the central bands for 7r + 7r~ 
recoiling against a selected p° or (j> and for K + K~ re- 
coiling against a p° are shown in Fig.0] For the /?°7r + 7r~ 
case we fit the ir + ir~ mass projection to the sum of a p° 
component, an /a(1270) component, and a p + p~ back- 
ground component. The p° is represented by the product 
of a P-wave relativistic Breit-Wigner with its width set to 
the Particle Data Group (PDG) value, a phase space 
term, and a factor 1 /rn^ w due to production via a virtual 



photon. The /2(1270) is represented by a D-wave rela- 
tivistic Breit-Wigner with its mean and width set to the 
PDG values. The p + p~ background shape is obtained 
from a sample of the related channel e + e~ — > p°p + p~ 
isolated by requiring two oppositely charged tracks iden- 
tified as muons. For the </>7r + 7r~ case, we use the same 
background parameterization in terms of fi and p + p~ , 
but refit for their normalizations. For the p°K + K~ case, 
we fit the K + K~ mass projection to the sum of a Breit- 
Wigner with mean and width fixed to their PDG values 
for the 4> signal, and a threshold function (q 3 )/ (1 + q 3 R), 
where q is kaon momentum in the cf> rest-frame and R is a 
shape parameter, for background. Assuming the masses 
of the two pairs to be uncorrelated and excluding the p° 
and <f> signal contributions, the fitted functions are inte- 
grated to obtain the tile fractions , ff, and /f . 

The extracted p°p° and 4>p° yields in the signal box 
are 1243 ± 43 and 147 ± 13 events, to be compared with 
total of 1508 ir+n-ir+TT- (~ 18% background) and 163 
K + K~it + it~ (~ 10% background) events in the signal 
box, respectively. 

To investigate the possibility of p°p° and <f>p° produc- 
tion in Y(4S) decay, we examine the data recorded at and 
below the T(4S) resonance separately. The yields below 
the T(45) resonance are 104 ± 14 for p°p° and 14 ± 4 for 
(j>p°, consistent with the expected values of 112 ± 4 and 
13± 1 obtained by scaling the on-peak yields of 1138±42 
and 135 ± 13 by the relative integrated luminosities. 

To investigate the production mechanism, we examine 
the production angle 9*, defined as the angle between 
the p® {4>) direction and the e~ beam direction in the 
CM frame. To measure the angular distributions, we 
subdivide the data into bins of 0* , and repeat the above 
fit, with linear background slopes s p o and s p o (s p o and 

s^) fixed to the values from the overall fit. The | cos6>*| 
distributions after MC efficiency correction are shown in 
Fig. The measurements are restricted to the fiducial 
region |cos6>*| < 0.8, as the efficiency drops rapidly be- 
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FIG. 4: Mass distribution for a) tt + tt~ pairs in p°iy + tv~ 
events, b) ir + n~ pairs in (pn + TT~ events, and c) K + K~ pairs 
in p°K + K~ events. The solid curves are the total fit. For 
the 7T + 7r - cases, the dotted curve is the p + p~ component, 
while the sum of /2(1270) and p + p~ contributions are shown 
as dashed. For the K + K~ case, the dashed curve represents 
the threshold function. 



yond 0.8. These forward peaking cos 9* distributions are 
consistent with the TVPA expectation which we find can 
be approximated by: 



da 



1 



cos 



g?cos 9* 



(3) 



in the fiducial region. The TVPA hypothesis gives a 
X 2 /dof (degrees of freedom) of 11.8/7 (p°p°) and 3.5/3 
{<j>p°). The fits disfavor 1 + cos 2 0*, giving a x 2 /dof of 
112/7 for p°p° and 6.3/3 for </>p°. 

Other observables are the <fi (p°) decay helicity angles 
9h, defined as the angle, measured in the <j> (p°) rest 
frame, between the positively charged kaon or pion and 
the flight direction of the <fi or p° in the CM frame. The 
efficiency-corrected distribution of cos9h , obtained using 
the procedure outline above for 9*, is shown for the p° 
and <j> candidates in Fig. |S| The solid lines in Fig. 
are normalized sin 2 6h distributions which give x 2 / dof 
of 19.3/9 0° from p°p°), 16.4/9 (</> from (f>p ), and 3.1/9 
(p° from 4>p°). The sin 2 9h distributions indicate that <f> 
and p° are transversely polarized as expected for TVPA. 
The dihedral angles, the angles between the decay planes 
of the two vector mesons measured in the CM frame, are 
consistent with a flat distribution with x 2 /dof of 7.0/9 
(p°p°) and 10.9/9 (<f>p°). 

The combined hardware and software trigger efficien- 
cies for signal events in the fiducial region are 99.9% for 



FIG. 5: Production angle distributions, after correction for 
efficiency, for a) p°p° and b) <j>p° '. The solid and dashed lines 
are the normalized i+ cos 9 an( } \ + cos 2 Q* distributions, re- 

1 — cos* 2 9* 

spectively. 



p°p° and 91.3% for <pp . The lower efficiency for cf>p° is 
due to an event shape cut in the software trigger. For the 
determination of signal cross sections, the MC cos#* and 
cos 9h distributions for <j> and p° are re- weighted to repro- 
duce the expectation from TVPA. The signal efficiencies 
in the fiducial region of |cos6**| <0.8 for p°p° and (j>p° 
are estimated to be 26.7% and 23.2%, respectively, in- 
cluding corrections to MC simulations of PID, tracking, 
hardware and software trigger efficiencies. Initial state 
photon radiation is included in the MC simulation. 

Systematic uncertainties due to PID and tracking ef- 
ficiency are estimated based on measurements from con- 
trol data samples. The related systematic uncertainties 
on lepton vetoes are estimated by the difference from not 
applying the e and p vetoes on pions. The systematic un- 
certainty from background subtraction is estimated by 
varying assumptions about background shapes. We in- 
vestigated possible feed-down background from related 
modes with an extra ir° using various extrapolations from 
the four-particle mass sidebands. We assume that the 
final states are fully transversely polarized. The system- 
atic uncertainties are summarized in Table [J 

Taking the branching fraction of <\> -^K + K~ as 49.1% 
and p° -^tt + it~ as 100% 3|, and signal mass regions of 
0.5 < m p0 < 1.1 GeV/c 2 and 1.008 < m < 1.035 GeV/c 2 , 
we obtain the following results for the TVPA cross sec- 
tions within | cos 9* \ <0.8 near s/s — 10.58 GeV: 



Cfid(e 



P°P° 



20.7 ± 0.7(stat) ± 2.7(syst) fb 
<°) = 5.7±0.5(stat)±0.8(syst) fb. 



The measured cross sections are in good agreement with 
the calculation from a vector-dominance two-photon ex- 
change model 0. 

In summary, we have observed exclusive production 
of C = +1 final states in e + e~ interactions. The mea- 
sured C parity configuration, the signal yields in data 
samples on the T(AS) resonance and below, and the pro- 
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FIG. 6: Decay helicity angle distributions for a) p° from p°p°, 
b) <fi from (ftp , c) p° from 4>P°- The solid lines are the nor- 
malized sin 2 9h distributions. 



TABLE I: Systematic uncertainties on the cross sections for 
e+e" ^pV/#°. 





pV 


# U 


Particle Identification 


9.6% 


10.4% 


Background subtraction 


7.0% 


7.0% 


Tracking efficiency 


5.0% 


5.0% 


p°p TY ,(/>p ( \ background 


1.6% 


2.7% 


Luminosity 


1.2% 


1.2% 


Total 


13.0% 


14.0% 



duction angle distributions support the conclusion that 
the production mechanism is two-virtual-photon annihi- 
lation. The Standard Model predictions of the anomalous 
magnetic moment of the muon and the QED coupling 
rely on the measurements of low-energy e + e~ hadronic 
cross sections, which are assumed to be entirely due to 
single-photon exchange. We have estimated the effect 
due to the TVPA processes we have measured, and find 
it to be small compared with the current precision 
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